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Abstract 

We consider a control problem for longitudinal vibrations of a non- 
homogeneous bar with clamped ends. The vibrations of the bar are 
controlled by an external force which is distributed along the bar. For 
' the minimization problem of mean square deviation of the bar we prove 

. that the optimal control has an infinite number of switchings in a finite 

' time interval, i.e., the optimal control is the chattering control. 

(N ■ 
O 

; 1. Introduction 

Consider small longitudinal vibrations of a nonhomogeneous bar of length /. 
^ \ The longitudinal displacement at a typical point x is denoted y{t,x) where 

if: is the time. Let g{x,t) be a density of external longitudinal force at the 
instant of time t at the point x. Suppose that 

g{t,x) = u{t)f{x) 

where the force profile function f{x) is assumed to be given, u (t) is the 
control function. We assume that 



-l<u{t)<l (1.1) 



*This research was supported by the Russian Foundation for Basic Research, 
project No. 11-01-00986-a . 
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The equation of longitudinal vibrations of the bar can be written as 

p (x) vttit, x) - {k (x) y^{t, x))^ = u{t)f{x) (1.2) 

Here 

p(x) = p (x) S (x) , k (x) = E (x) S (x) 

where p (x) is the density of bar, 5* (x) is the cross-sectional area, E (x) is 
the Young's modulus at x, see, for example, [11 [21 [3]. 
We assume that the ends of the bar are clamped: 

l/Uo=2/U^ = 0, t>0 (1.3) 
and the initial position and velocity are fixed: 

y\t=o = yoix), xe[o,i] (i.4) 



yt\t=o = ^ ^ [0'^] (1-5) 

We suppose that the coefficient functions k, p are smooth enough and 

Vxe[0,/] A;(x)>A;o>0, p{x)>po>0 (1.6) 

We consider an optimal control problem: to find such a control function 
u (t) that minimize the following functional 

POO rl 

/ / p (x) {t, x) dxdt ini (1.7) 
Jo Jo 

under ffT:2|) -f fL6|) . 

The problems of longitudinal vibrations of a bar were considered in [H [21 
[3]. In [H l2j the dynamics of the longitudinal vibrations of a bar subjected 
to viscous boundary conditions was studied. The optimal boundary control 
problem for longitudinal vibrations of a bar was considered in [3J . By using a 
maximum principle the optimal control was expressed in terms of an adjoint 
variable. 

In this paper for the problem of controlling the longitudinal vibrations of 
a bar f ll.ip -f lTTTl) we construct a solution y {t, x) in the form 

oo 

yit,x) = J2s,it)hjix) (1.8) 
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where {hj are eigenfunctions of the Sturm-Liouville problem, {sj (^)}°^ 

are corresponding Fourier coefficients. To find Fourier coefficients we consider 
an optimal control problem in the space For the control problem in we 
show that the optimal solutions contain singular trajectories and chattering 
trajectories. A trajectory is called a chattering trajectory if it has an infinite 
number of a control switchings on a finite time interval. By similar method 
we studied the optimal control problem for a rotating uniform Timoshenko 
beam [U [5]. But for the Timoshenko beam similar results hold only for a 
dense set of initial conditions in the space 

2. Optimal control problem in I2 

Define an operator L in ([0, /]) by 

Lh = (kh^)^ 

Consider the following Sturm-Liouville eigenvalue problem with Dirichlet 
boundary conditions: 

Lh + Xp{x)h = 0, xe{0,l) (2.1) 

h{0) = 0, h{l)=0 (2.2) 

Here the functions k (x) and p (x) satisfy (II. 6p . It is known (see |6j) that 
the problem (I2.1l) - (l2.2p has an infinite sequence of eigenvalues {Aj}°^^, which 
are simple and positive: 

< Ai < A2 < . . . , Aj — )■ 00, j ^ 00 

To each eigenvalue Xj corresponds a single eigenfunction hj, and the sequence 
of eigenfunctions {hj{x)}'^^ forms an ortonormal basis of L2 {{0,1) ;p) with 
the inner product {z, w)^ = SqP {x) z {x) w {x) dx. If k{x) ,p {x) are smooth 
enough (for example, k' , p E ([0,/])) and the condition (11. 6p holds, then 
the eigenvalues Xj admit the asymptotic form [3 [HI [9]: 

^ ~ TT^ (j^ ^p{x) /k {x)di\ , J ^ 00 (2.3) 
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Assume that y{t, ■) G L2 ((0, /) For any t > we expand the solution 
y {t, x) of fll.2p in the basis {hj{x)}'^y 

y{t,x) = (2-4) 
Sj{t) = / p{x)y{t,x)hj{x)dx = {y,hj)^ 

J 

Using ([HD-dES]) we get: 

"Sj (0) = / P (x) yo (x) hj (x) dx = (yo, hj) 



■Sj (0) = / P (x) yi (x) hj (x) dx = {yi, hj)^ 
We multiply the equation (ll.2p by hj and integrate it in x: 



{pytt + Ly) hjdx = ufhjdx 



d^ d^ 

^ (y^ hj)p + {Ly, hj) = u if, hj) — {y, hj)^ + {y, Lhj) = u (/, hj) 

or 

d^ 

(y, hj)p + \j {y, hj)^ = u (/, hj) 
Thus the function Sj (t) satisfies the following equation: 
Sj{t) + \s,{t) = C,u{t), J = 1,2,... 

where 

Cj = {f,h,)= f f{x)hj{x)dx (2.5) 
We substitute (12.41) into fll.7p . Using Parseval's equality we get: 

roo rl POO °° 

/ / p{x)y^ {t,x)dxdt = / Vs.(t)o?t (2.6) 
Jo Jo Jo j^^ ■' 

Denote aj = (?/o,^j)p, /3j = ivii^j)^. We reduce the problem ( 11.2p -( fL7ll 
to the following one: 

POO 

/ y^s]{t)dt^M (2.7) 
Jo ■' 



L. Manita 



Optimal control for longitudinal vibrations 



Sj{t) + \,Sjit) = CMt), J = 1,2,... (2.8) 

Sj{0) = aj, Sj{0) = (3j (2.9) 

- 1 < u{t) < 1 (2.10) 
We shall assume everywhere below that 

Cj^O for all j = 1,2,... (2.11) 

Remark. Assumption (12. lip is very essential for the problem (I2.7p - (l2.10p . 
Indeed, let Cj^ = for some jo- Then jo-th equation in (12. 8p takes the form 

Sjoit) + KSj,{t) =0 

Hence, if |a;jo| + ^ then the corresponding solution sj^ (t) does not 
vanish as t — )■ oo. Therefore the integral (12. 7p is equal to +oo and the 
optimization problem (l2.7p -f l?!TU]) has not any sense. 

Assume that 

?/o, z/iGi^2((o,0;p), /eL2(o,0 (2.12) 

Following |5] we denote 

= \Ai ' ^ji^) = Sj{t) /ujj, Cj = Cj/ujj, aj = aj, bj = (3j/ujj 
Then the problem (I2.7l) - (l2.10p takes the form 

I]s](t)dt ^ inf (2.13) 

Sj = UjTj, fj = -UjSj + CjU (2.14) 

s,(0) = a„ Tj{0) = b,, J = 1,2,... (2.15) 

- 1 < u{t) < 1 (2.16) 



Denote 



sit) = (si(t),S2(t),...), r(t) = (ri(t),r2(t),...) 
a = (ai,a2, . . .) , 6 = (61, 62, • • •) , c = (ci, C2, . . .) 
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Consider the standart Hilbert space l2- 

^2 = |w = (wi,W2, • • •) : w„ G R, ^^w^ < ooj 

with inner product {v, w) = VnWn- Using assumption f l2.12p we get that 
a, b, c E l2- 

The existence and uniqueness of a solution (s (t) ,t {{)) to problem (I2.13p - 
fl2.16p in the space I2 x I2 were proved in [S] for any initial data from an open 
neighborhood of the origin (s = 0, r = 0). We apply a formal generalization 
of the Pontryagin maximum principle to the problem fl2.13p - fl2.16p . Denote 
by ipi = {il^ii,ipi2, ■ ■ ■) {i = 1,2) adjoint variables. Define the Pontryagin 
function 

00 
i=i 

where 

^^0 (^1, 'ip2, S,t) =Y^ [i^lj^n - i^2j^jSj Hi {tpi,1p2, S,t) = Y, i^2jCj 

i=i ^ ^ ^ 3=1 

For breavity we denote z = {ipi,ilj2, s, r). In the space I2 x I2 x h >^ h let us 
consider the Hamiltonian system 

^ij = i^2ji^j + Sj, Sj = UjTj J = 1 2 (2 17) 

where u*{t) satisfies the following maximum condition: 

u*{t) = argmax H {z (t) , u) = argmax {uHi {z (t))) (2-18) 

Here we use notation: a* = argmax 5f(a) iff g{a*) = ma.xg{a). 

It was proved [5j that the Pontryagin maximum principle is the necessary 
and sufficient condition of optimality for the problem fl2.13p -f l2.16p . 

If Hi[z{t)) 7^ along the trajectory the optimal control is uniquely 
determined as a function of time from the maximum condition fl2.18p : 

(00 
Y^p2j{t)Cj 
3 = 1 
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Suppose that there exists an interval (^1,^2) such that 

Hiiz{t)) = 0, Vte(ti,t2) 

To find an optimal control u (t) in this case we will differentiate the identity 
Hi {z (t)) = by virtue of the system (I2.17P until the control u with a non- 
zero coefficient occurs in the resulting expression with a non-zero coefficient. 



dt 

d^ 
dt^ 



J2T7t 



HAz) 



HAz) 



d 
di 

d_ 
dt 



1213 L?/''''' 



00 



Ui] — U > CaUJ 



2 2 



(2.19) 



|2.17t j=i 

Assume that all series in (I2.19P are convergent in t^. Denote 

00 00 

H2 (-2) = - E ^^li'^i' ^3 (^) = - E ^3^3 (^sjtUj + Sj) 
j=i i=i 

00 

^4 (^) = - E Cj'^j (-V'ljWj + 'Tj) 

From fl2:T9|) it follows that 

Hi {z (t)) = H2 {z (t)) = H^ [z (t)) = {z (t)) = 0, t G (ti, h) 
We say a solution of the f l2.17p -( l2.18p is singular if it belongs to the surface 
T. = {z: Hi{z) = H2iz) = Hsiz) = H^iz) = 0} (2.20) 



A singular control (t) is determited from the equation 
Using f l2.19p we obtain 



d4 



I2T71 



u 



°(t) = ^c,a;|(^2,c^,+2s,)/Ec: 



2 2 

'3^j 



HAz) = 0. 



(2.21) 
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Note that the origin [ipi =0, ip2 = 0, s = 0, r = 0) is the singular trajectory 
and the corresponding singular control {t) equals 0. 

It was proved [5j that in a certain neighborhood of the origin the structure 
of the optimal solutions is the following one: for the finite time the optimal 
nonsingular trajectory enters the singular surface with infinite numbers of 
control switchings, after that the optimal trajectory remains on the singular 
surface and attains the origin for the infinite time. Namely, the following 
theorem holds. 

Theorem 1 [5j. Let cj ^ Vj and (ciwf , C2C<;2) . . .) € li. Assume 
that there exist positive constants b and K such that 

Iwj+il - > 5, \ujj\<K-2, j = l,2, ... 

Then there exists an open neighborhood of the origin in the space (s, r) 
such that the following statements hold for all initial data (a, 6) from this 
neighborhood. 

(i) The problem fl2.13p -( !2.16p has a unique optimal solution. 

(ii) In the space z = {ipi, ip2, s,t) there exists the singular surface S of 
codimension 4 given by the equations 

00 00 

j2 'ip2jCj = 0, Cji^ij^j = 

00 00 

^3^0 i'p23^3 + S3) =0' H Cj^j i-^13^3 + ^3) = 

i=l 3=1 

which is filled in by the singular extremals of the problem (I2.13p -( l2.16p . 
The control on singular extremals are defined by (I2.2ip . 

(iii) For all initial data not belonging to the projection of the singular surface 
S on the space (s, r), the optimal trajectories arrive at S in finite time 
with countable many control switchings, i.e., the optimal trajectories 
are chattering trajectories. 
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3. Optimal solution for controlling vibrations 

Let {s* {t) ,u* (t)) be an optimal solution of (12.13p -( j2.16p . Consider 



where {hj{x)}°°^^ are eigenfunctions of the Sturm-Liouville problem (12. ip - 



Series (13. ip formally satisfies equation (II. 2p . boundary conditions (II. Sp 
and initial conditions ( ll.4p - (ll.5p . We will show that this series gives a weak 
solution of problem (ll.2p - (ll.5p . 

Denote Qt = {0,1) x (0,T), where T > 0. Consider the Sobolev space 
H'' (Qt) = W2 (Qt) , k > 0. The space (Qt) consists of all functions 
V (z L2 {Qt) whose generalized derivatives up to order k exist and belong to 
L2 {Qt)- The space {Qt) can be defined as a completion of {Qt) with 
respect to the norm of the space {Qt)- 

Let g = uf e L2{Qt) , |/i 6X2(0, 

Definition 1 . We say a function y G {Qt) is a weak solution of 
the (ir2D-([L5]) if 

(i) yL=o = y\x=i = 0) ^ > 0) y\t=o = yo{x), x e [o, /]; 



for each v e {Qt)- v \x=o = v \x=i = 0, v \t=T = 0. 

Definition 2 . We say y G {Qt) is a almost everywhere solution of the 
problem (ll.2p -( TT3]) provided y satisfies equation ( II. 2p in Qt for almost all 
{t,x) G Qt and y satisfies (ll.3p - (ll.5p . 

The following theorem is the main result for the problem (ll.2p - (ll.7p . 

Theorem 2. Let ?/o £ Hq{0,1), yi G Hq{0,1), k{x),p{x) are smooth 
enough (for example, k, p E C'^{[0,1])), p{x) > Po > 0, k{x) > /cq > 0. 



Assume that / G C^[0,/], 

/(0) = /(/)=0, /W(0) = /» (0=0, 2 = 1,2,3 (3.2) 



00 



y* {t,x) = Y.s*{t) h,{x) 



(3.1) 



(12:2]). 
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and condition (12.1 ip holds. Then there exist positive constants qi and q2 
such that if 

\\yo\\L2miy,p) < II^iIIl2((o,o;p) < ^2 

then 

(i) the problem fll.2p - (ll.7p has a unique optimal solution y* {t,x); 

(ii) y* G (Qt) for all T > 0; 

(iii) an optimal solution y* (t, x) has an infinite number of control switchings 
in a finite time interval. 

Proof. Here we use notations introduced in Section 2. Since the functions 
p, k, f satisfy conditions of Theorem[2]it follows [S] that Cj ~ as j — )■ oo, 
where ^ 

Cj = (/, hj) = f (x) hj (x) dx 

J 

Then we have 

E (ci^J) = E (Cj^f ) = E ^1 < ^ ^ (ci^^ cs^t • • •) e /2 

i=i i=i i=i 

The property (12.31) of the eigenvalues {Aj}°^-^ of the problem fl2.ip - fl2.2p 
imply that there exist positive constants 5 and B such that 

IcUj+il — l^jl > 5, l^jl < i?j 

Now we may apply Theorem [T]to the problem fl2.13p -f l2.16p . We get that the 
optimal control u* (t) for the problem fl2.13p -f l2.16p has an infinite number 
of switchings in the finite time interval. 

Since the functions /, a, /3 satisfy the conditions of Theorem 12] it follows 
(see [TU| ITT] ) that the function y*{t,x) defined by (13. ip is a unique weak 
solution of the problem flL2|) -f ir5|) and y* G H'^ (Qt)- Hence [lOj y*{t,x) is 
the solution almost everywhere of the problem fll.2p - fll.5p . Thus the function 
y*{t,x) satisfies (11. 2p for almost all {t,x) G (0,/) x (0, +00), boundary and 
initial conditions fll.3p -f lT3]) . 

Since the function s*{t) = (sj(t), S2{t), . . .) minimizes the functional fl2.7l) 
and the identity (12. 6p holds then the function y* {t, x) minimizes the functio- 
nal (11. 7p . Thus y* {t, x) is a solution of the problem fll.ip -f lL7|) . 
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4. Conclusion 

We considered the optimal control problem of longitudinal vibrations of a 
nonhomogeneous bar with clamped ends. We proved that the optimal tra- 
jectories contain singular part and nonsingular one with accumulation of 
control swithings. 
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